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Despite our increasing topological knowledge on regulation networks in model bacteria, it is largely 
unknown which of the many co-occurring regulatory events actually control metabolic function and 
the distribution of intracellular fluxes. Here, we unravel condition-dependent transcriptional 
control of Escherichia coli metabolism by large-scale 13 C-flux analysis in 91 transcriptional regulator 
mutants on glucose and galactose. In contrast to the canonical respiro-fermentative glucose 
metabolism, fully respiratory galactose metabolism depends exclusively on the phosphoenol- 
pyruvate (PEP)-glyoxylate cycle. While 2/3 of the regulators directly or indirectly affected absolute 
flux rates, the partitioning between different pathways remained largely stable with transcriptional 
control focusing primarily on the acetyl-CoA branch point. Flux distribution control was achieved 
by nine transcription factors on glucose, including ArcA, Fur, PdhR, IHF A and IHF B, but was 
exclusively mediated by the cAMP-dependent Crp regulation of the PEP-glyoxylate cycle flux on 
galactose. Five further transcription factors affected this flux only indirectly through cAMP and Crp 
by increasing the galactose uptake rate. Thus, E. coli actively limits its galactose catabolism at the 
expense of otherwise possible faster growth. 
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Introduction 

To survive environmental vicissitudes, microbes evolved 
regulatory systems that sense external (and internal) condi- 
tions to mount appropriate cellular responses (Janga et al, 
2007; Kotte et al 2010). A large fraction of these functional 
responses concern metabolic adaption to changing nutritional 
conditions. While activity of the metabolic enzyme is subject 
to multiple levels of regulation, transcriptional regulation is 
considered a main mechanism for such metabolic adaptations 
in bacteria (Gama-Castro et al, 2008), but there are indications 
that its relevance might have been overemphasized (Heine- 
mann and Sauer, 2010). Understanding the regulatory pro- 
cesses that govern cellular metabolism has become a key focus 
in microbial systems biology. Key elements of such systems 
biology approaches are genome-scale models of metabolic 
stoichiometry (Feist and Palsson, 2008) and the detailed topology 
of transcriptional regulatory networks that describe all known 



interactions between transcription factors and their target genes 
(Gama-Castro et al, 2008). Unfortunately, such topological 
networks offer only a static description of potentially occurring 
interactions, of which only a subset will be active at any given 
point in time and conditions (Harbison et al, 2004; Luscombe 
et al, 2004; Balazsi et al, 2005; Moxley et al, 2009). 

This condition dependence of transcriptional regulation is 
typically assessed by measurements of mRNA abundance 
changes. Typically unanswered remains the question which of 
all the active regulation processes actually control function, in 
particular metabolic function, and to which extent? To assess 
metabolic function in the network context, intracellular 
reaction rates (i.e., fluxes) must be quantified by methods of 
13 C-based metabolic flux analysis (Sauer, 2006) . By comparing 
the network's flux response with regulatory events, one can 
identify actual control of function. In some specific cases this 
has been done, e.g., by demonstrating that the redox regulator 
ArcA controls tricarboxylic acid (TCA) cycle fluxes during 
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aerobic batch growth on glucose (Perrenoud and Sauer, 2005), 
but exerts no apparent flux control during glucose-limited 
growth of Escherichia coli (Nanchen et al, 2008) . The cAMP 
receptor protein Crp on the other hand specifically controls 
phosphoenol-pyruvate (PEP) -gly oxylate cycle fluxes during 
glucose-limited growth of E. coli, but not during excess batch 
growth (Perrenoud and Sauer, 2005; Nanchen et al, 2008). 
So far, only one larger-scale flux study has recently been 
concluded, revealing several condition-specific networks of 
transcriptional regulation that control metabolic function in 
yeast (Fendt et al, 2010) . Overall, these results clearly demon- 
strate the strong condition dependence of metabolic control with 
typically sparse networks of active transcriptional control that 
control the distribution of fluxes between different pathways. 

Focusing on central carbon metabolism of E. coli, we aim 
here to systematically identify transcriptional regulators that 
control the distribution of metabolic fluxes during aerobic 
growth on hexoses. For this purpose we selected 81 transcrip- 
tion factors, including all known to directly or indirectly 
control central metabolic enzymes, and 10 sigma- and anti- 
sigma factors. Of the 81 transcription factors, 41 have one or 
more direct target gene that is involved in central carbon 
metabolism (Supplementary Table 1). To assess the condition 
dependence of transcriptional control of flux, we selected 
glucose and galactose as two substrates that are highly similar, 
yet lead to distinct growth rates (Soupene et al, 2003), overall 
metabolic rates (De Anda et al, 2006; Samir El et al, 2009) and 
levels of catabolite repression (Hogema et al, 1998; Betten- 
brock et al, 2007) . Different from the recent yeast flux study 
(Fendt et al, 2010), we quantify here not only flux ratios 
but also absolute values for intracellular fluxes by large-scale 
13 C-constrained flux analysis (Zamboni et al, 2009). From 
these data, we show here that E. coli employs substantially 
different modes of metabolic operation and pathway usage on 
both substrates, and identify the controlling transcription 
factors and their target pathways. 

Results 

Two distinct modes of hexose catabolism on 
galactose and glucose 

To quantify physiological and intracellular flux difference of 
aerobic glucose and galactose catabolism, E. coli wild type was 
grown in separate shake flask experiments with 3 g l -1 of either 
the [1- 13 C] hexose isotope isomere or a mixture of 20 % [U- 13 C] 
with 80% natural abundance hexose isotope isomers. 
Compared with the glucose-grown culture, galactose-grown 
cultures grew substantially slower with an about fourfold lower 
metabolic rate and almost absent overflow to acetate (Table I). 
To elucidate how these different macroscopic rates were 
reflected by the intracellular distribution of fluxes, we first 
determined ratios of converging fluxes from GC-MS- detected 
mass isotope partitioning in proteinogenic amino acids (Table I) 
(Zamboni et al, 2009). These ratios indicated major alterations 
in the flux connecting oxaloacetate and PEP as well as in the 
contribution of the glyoxylate cycle to the oxaloacetate pool. 

For a better resolution of these flux differences, we 
estimated absolute intracellular fluxes using the physiological 
data and the flux ratios of Table I as constraints (Zamboni et al, 



Table I Physiological parameters and ratios of converging fluxes during 
wild-type growth on glucose or galactose 

Glucose Galactose 



Wild-type physiology* 
Growth rate (h~ Y ) 

Hexose uptake (mmolgCDW -1 h _1 ) 
Acetate secretion (mmolgCDW -1 h -1 ) 

Wild-type flux ratios 0 
EMD, ED and PP pathways 
Serine through EMP pathway 
Pyruvate through ED pathway 
PEP from PP pathway (upper bound) 



0.61 ±0.01 0.18 ±0.01 
8.26 ±0.50 2.00 ±0.33 
4.89 ± 1.52 0.11 ±0.16 



0.77 ±0.01 0.79 ±0.01 
0.07 ±0.02 0.04 ±0.05 
0.20 ±0.07 0.35 ±0.07 



Gluconeogenesis, TCA and glyoxylate shunt 

Oxaloacetate from phoshoenolpyruvate 0.68 ± 0.01 0.20 ± 0.03 

Oxaloacetate from glyoxylate shunt NA 0.72 ± 0.20 

Pyruvate from malate (upper bound) 0.04 ± 0.02 0.02 ± 0.03 

Pyruvate from malate (lower bound) 0.01 ± 0.00 0.02 ± 0.03 

Phoshoenolpyruvate from oxaloacetate 0.02 ± 0.00 0.20 ± 0.03 



CI metabolism 

Serine from glycine 
Glycine from serine 



0.31 ±0.02 
0.98 ±0.02 



0.46 ±0.02 
0.99 ±0.02 



a Values and standard deviations were obtained from at least three biological 
replicates. 

b Values are flux ratio values ± 95 % confidence intervals. Two experiments led to 
identical results. 



2009; Figure 1), and confirmed the results with whole 
isotopologue modeling (Kleijn et al, 2005; Supplementary 
Tables 2 and 3). Consistent with earlier data (Fischer and 
Sauer, 2003b; Perrenoud and Sauer, 2005), we found the 
classical glucose flux distribution with ~30% of the high 
hexose influx channeled into the pentose-phosphate (PP) 
pathway and acetate overflow that exceeded the TCA cycle flux 
about fourfold. The relative distribution of the much lower 
galactose influx was similar to glucose metabolism in the 
upper part of metabolism. In the lower part, however, the 
distribution was entirely different with a substantial glyoxylate 
shunt flux and an interrupted TCA cycle between 2-oxogluta- 
rate and succinate. Conjointly with a higher flux through PEP 
carboxykinase, galactose metabolism, hence, relies almost 
exclusively on the so-called PEP-glyoxylate cycle for respira- 
tory energy generation (Fischer and Sauer, 2003b). Specifi- 
cally, two PEP molecules are converted to two acetyl-CoA, 
one of which enters the glyoxalate shunt reactions to form 
succinate and the other is fused with glyoxlate to form malate. 
Both organic acids are converted to oxaloacetate, one of 
which enters the next round of the cycle and the other is 
decarboxylated to form PEP. With a stoichiometry of 2 PEP -> 
3C0 2 + 4NADH + UQH 2 + ATP + PEP, this PEP-glyoxylate 
cycle hence bypasses the cyclic operation of the TCA cycle by 
the joint activity of PEP carboxykinase and glyoxylate shunt, 
such that their normally considered functions in gluconeogen- 
esis and anaplerosis, respectively, are recruited for complete 
combustion of hexoses to C0 2 (Fischer and Sauer, 2003b) . 

Intracellular flux partitioning in 91 transcription 
factor mutants on glucose and galactose 

Having established two different modes of hexose catabo- 
lism — rapid respiro-fermentative metabolism on glucose and 



2 Molecular Systems Biology 2011 



© 2011 EMB0 and Macmillan Publishers Limited 



Large-scale 13 C-fluxes in E. coli regulator KOs 
BRB Haverkorn van Rijsewijk et al 



B GALACTOSE 




TCA CYCLE & 
GLYOXYLATE 129 
SHUNT 



SUC-CoA^_\J&^ 2-OG 

Figure 1 Absolute metabolic fluxes in E. coli during aerobic growth on glucose (A) or galactose (B). Flux arrows are drawn in proportion to the substrate uptake rates 
for each condition. The numbers represent absolute flux values (mmolgCDW -1 h" 1 ). One of two replicate experiments is shown (Supplementary Tables 2 and 3). 
The presented fluxes are from one of two independent experiments and were obtained by 13 C-constrained flux analysis using the software FiatFlux (Zamboni etal, 
2005). They were independently confirmed by flux estimation with a whole isotopologue model (Kleijn et al, 2005; van Winden et al, 2005). Generally, the deviation 
between the two independent experiments was 1-5%, and whole isotopologue sensitivity analysis through addition of Gaussian noise confirmed accurate estimation for 
all major fluxes (Supplementary Tables 2 and 3). 



slower fully respiratory metabolism via the PEP-glyoxylate 
cycle — we next wanted to understand how these modes were 
transcriptionally controlled. For this purpose, a large-scale flux 
analysis was performed with transcription factor knockout 
mutants under both growth conditions. Specifically, we chose 
81 mutants of transcription factors that are known to directly 
or indirectly control enzymes of central metabolism (Keseler 
et al, 2005) and 10 sigma- and anti-sigma factors that control 
expression globally (Supplementary Table 1). Fully aerobic 
deep-well microtiter plate cultures (Duetz et al, 2000) were 
grown on either the 100% [1- 13 C] -labeled hexose isotope 
isomere or a mixture of 20% (wt/wt) [U- 13 C] and 80% natural 
abundance hexose (Fischer and Sauer, 2003a). Physiology and 
flux distribution of wild-type deep-well cultures were highly 
similar to the above shake flask cultures (Supplementary 
Tables 2 and 3). For each deletion mutant, both independently 
generated clones in the library were analyzed separately 
and yielded similar results, except for the Crp and Cra 
mutants. PCR verification revealed that both genes were still 
present in one of the two clones, which were therefore 
discarded. 

On the basis of 13 C-determined flux ratios and extracellular 
rates, we estimated absolute intracellular fluxes for all mutants 
(Supplementary Tables 2 and 3; Zamboni et al, 2005). The 
average deviation of absolute fluxes through the 23 central 
metabolic reactions for all 91 mutants from the wild-type value 
was 19 and 49% on glucose and galactose, respectively 
(Figure 2). In contrast to absolute fluxes, the relative flux 
partitioning was rather invariant in the entire mutant set. 
Genetic perturbations hence primarily affected absolute fluxes 
and not the distribution of fluxes as was described before 
(Blank et al> 2005; Fischer and Sauer, 2005; Perrenoud and 
Sauer, 2005; Ishii et al, 2007; Nanchen et al 2008). 
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Figure 2 Average deviations of absolute fluxes and flux ratios in 91 mutants 
from the wild type. Deviations in absolute fluxes and flux ratios are %-differences 
compared with the wild-type values for each condition. 

Next we focussed on flux changes that occurred for key 
metabolic pathways in the 91 mutants (Figure 3). On glucose, 
transcription factor deletions altered many fluxes, including 
those of hexose uptake, glycolysis, pyruvate dehydrogenase, 
TCA cycle and acetate secretion (Figure 3A). On galactose, 
these flux changes were even larger and additionally included 
the glyoxylate shunt. These variations in absolute flux values 
may be caused by either a change of the overall flux magnitude 
(due to or resulting in increased carbon uptake) or changes in 
the distribution of flux within the network; hence, affecting 
specific flux magnitudes. To distinguish between these 
two possibilities, we also compared fluxes normalized to the 
uptake rate between wild type and the 91 mutants (Figure 3B) . 
Many transcription factors affected absolute uptake rates. 
About 60 of the 91 transcription factor mutants exhibited more 
than a 10 % deviation in uptake rates compared with the wild 
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Figure 3 Absolute (A) and substrate uptake normalized (B) flux changes of key metabolic pathways in the 91 mutants compared with the wild type during growth on 
glucose (A) and galactose (□). The dashed line indicates the wild-type reference fluxes. 



type under either condition, indicating the importance of 
transcriptional regulation in controlling metabolic flux. 
The distribution of key fluxes through glycolysis, PP and 
Entner-Doudoroff (ED) pathway in upper metabolism, how- 
ever, remained constant under both conditions (Figure 3B). 
Thus, flux variations in upper metabolism were caused by an 
altered magnitude of overall fluxes, but our data do not allow 
differentiating between a putative-specific transcriptional 
control of pathways in upper metabolism and indirect, 
growth-related effects. In key fluxes of lower metabolism 
(pyruvate dehydrogenase, anaplerosis, TCA cycle, glyoxylate 
cycle, acetate secretion), not only absolute fluxes but also their 
distribution changed (Figure 3B). TCA cycle, glyoxylate shunt 
and acetate secretion fluxes normalized to uptake rates were 
significantly different in several mutants, indicating that in 
lower metabolism also the distribution of flux is transcription- 
ally controlled. 



Identification of transcription factors that control 
the distribution of flux at key branch points 

To identify individual transcription factors that actually 
control the distribution of fluxes within the network, we 



focused our attention on key metabolic branch points. As 
relative fluxes in upper metabolism remained constant 
(Figure 3B), the almost perfect linear correlation between 
hexose uptake and glycolytic flux at the glucose-6-phosphate 
branch point in the various transcription factor mutants was 
expected (Figure 4A). The only exceptions were control 
experiments with disruptions of either glycolysis or PP 
pathway in the Pgi and Zwf mutants, respectively. Clearly, 
transcriptional regulation does not control the distribution 
of flux into the PP pathway during growth on glucose or 
galactose. Furthermore, flux partitioning at the glucose-6- 
phosphate branch point was rigid over a large range of 
absolute catabolic fluxes from 1.0 to 9.0mmolg _1 h _1 . This 
rigidity in the face of substantially different biomass yields in 
the different mutants (Supplementary Tables 2 and 3) 
demonstrates that, during hexose degradation, the PP pathway 
flux is not determined by the biosynthetic needs for the redox 
cofactor NADPH or pentoses, as was previously described for 
Bacillus subtilis (Fischer and Sauer, 2005) . 

As the flux distribution varied significantly in lower 
metabolism (Figure 3B), we focused on the acetyl-CoA branch 
point in more detail. For the glucose-grown mutants we plotted 
the incoming pyruvate dehydrogenase flux as a function of the 
outgoing flux into the TCA cycle (i.e., citrate synthase flux; 
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Figure 4 Influence of individual transcription factors on the flux distribution at key metabolic branch points for growth on glucose (A) and galactose (■). The four 
branch points are the glucose-6-phosphate node on both hexoses (A), the acetyl-CoA node on glucose (B) and on galactose (C), and the isocitrate node on galactose 
(D). Values for unlabeled mutants are the average of two independent experiments obtained by 13 C-constrained flux analysis (Supplementary Tables 2 and 3). For 
mutants with significantly altered flux distributions (labeled), a third independent experiment was conducted and the error bars represent standard deviations of these 
three independent experiments. Enzyme control mutants are given in italics. The wild-type values are highlighted by black circles. 



Figure 4B) . Flux partitioning at this branch point was indeed 
transcriptionally controlled, because nine mutants exhibited 
significantly altered flux partitioning. For the mutants Fur 
(Escolar etal, 1998; Zhang etal, 2005), PdhR (Quail and Guest, 
1995; Ogasawara et al, 2007), ArcA (Iuchi and Lin, 1988; 
Perrenoud and Sauer, 2005; Vemuri et al, 2006) and IHFA and 
B (Bongaerts et al, 1995; Green et al, 1997; Park et al, 1997; 
Cunningham and Guest, 1998) altered partitioning can be 
explained by the known targets of these transcription factors in 
the TCA cycle and/or the respiratory chain. The different flux 
partitioning effects of the two IHF mutants can be explained by 
different stabilities of the IHF homodimer B/B and A/A 
complexes that have different capacities to restore the wild- 
type phenotype of the normal heterodimer complex (Werner 
et al, 1994; Zulianello et al, 1994). Altered flux partitioning 
between respiration and acetate overflow in the four remain- 
ing transcription factor mutants GlpR (glycerol metabolism; 
Larson et al, 1987), GlcC (glycolate and glyoxylate metabo- 
lism; Pellicer et al, 1996), QseB (quorum sensing; Sperandio 
et al, 2002) and HdfR (flagella formation; Ko and Park, 2000) 
was not overly strong, but significant. Without known target 
genes in respiration, TCA cycle or overflow metabolism, these 



mutant phenotypes might be caused through indirect effects or 
through yet unknown target genes of these transcription 
factors. Overall, of all 36 transcription factors known to 
regulate enzymes of the diverging pathways at the acetyl-CoA 
branch point (i.e., TCA cycle, glyoxylate cycle, acetate 
secretion), only a minority actively regulates flux during 
steady-state growth on glucose. 

In contrast to glucose, a third outgoing flux occurred at 
the acetyl-CoA branch point during galactose metabolism, i.e., 
fusion of acetyl-CoA with glyoxylate in the malate synthase 
reaction of the PEP-glyoxylate cycle. To differentiate potential 
regulatory events in either of the three outgoing pathways at 
the acetyl-CoA branch point (i.e., acetate secretion, TCA cycle 
flux and PEP-glyoxylate cycle flux), we plotted the incoming 
pyruvate dehydrogenase flux as a function of the summed 
outgoing flux into TCA and PEP-glyoxylate cycle (Figure 4C) 
and the flux from acetyl-CoA to isocitrate as a function of the 
PEP-glyoxylate cycle flux (Figure 4D). Several mutants 
exhibited an increased metabolic overflow with reduced 
overall flux into the TCA and PEP-glyoxylate cycle on galactose 
(Cra, IHFA, IHF B, Mlc, NagC). Deletion of these transcription 
factors also abolished or reduced the PEP-glyoxylate cycle, 
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redirecting flux through the TCA cycle instead (Figure 4D). 
While not affecting overflow metabolism, also the Crp mutant 
abolished the PEP-glyoxylate cycle. This Crp-dependent 
control of the PEP-glyoxylate cycle flux was previously 
described for slowly growing glucose-limited chemostat 
cultures, presumably triggered by high cAMP concentrations 
that allosterically activate Crp (Nanchen et al, 2008) . 



Increased hexose uptake rate results in 
altered cAMP-Crp-dependent control of 
the PEP-glyoxylate cycle 

In contrast to Crp, the activity of the other five transcription 
factors is not modulated by cAMP; hence, their influence on 
the PEP-glyoxylate cycle must be initiated by different means. 
As both Cra (Saier and Ramseier, 1996; Crasnier-Mednansky 
et al 1997) and Mlc (Decker et al, 1998; Kimata et al 1998; 
Plumbridge, 2002) are regulators of the glucose phospho- 
transferase system (PTS Glc ) that represses the uptake of non- 
PTS sugars such as galactose (Misko et al 1987) and as NagC is 
a known direct repressor of the galactose transporter GalP 
(Soupene et al 2003; Samir El et al 2009), we wondered 
whether there was a general relationship between galactose 
uptake and the PEP-glyoxylate cycle flux. Indeed, when 
plotting the fraction of acetyl-Coenzyme A that enters the 
PEP-glyoxylate cycle at the acetyl-Coenzyme A branch point 
against the galactose uptake rate for all mutants, the strongly 
reduced fraction of carbon entering the PEP-glyoxylate cycle in 
the Cra, IHF A, IHF B, Mlc and NagC mutants correlated with 
much higher galactose uptake rates than in the wild type 
(Figure 5 A) . As Mlc is actively sequestered to the membrane by 
the glucose PTS Glc (Plumbridge, 2002), we hypothesized that 
the PTS Glc might have an active role in the regulation of 
galactose uptake. We therefore grew the PTS Glc -component 
IIA glc mutant Crr on galactose. Consistent with the hypothesis 
of a direct or indirect role of the PTS Glc in galactose uptake, the 
galactose uptake rate of the Crr mutant was indeed signifi- 
cantly increased (Figure 5A). Additionally, the PEP-glyoxylate 
cycle was repressed in the Crr mutant, as was also seen for all 
other mutant with increased galactose uptake. 

To verify causality of the observed correlation between 
hexose uptake rate and PEP-glyoxylate cycle usage, we used an 
environmental strategy to decrease hexose uptake. Specifi- 
cally, we performed glucose- and galactose-limited continuous 
culture experiments at different dilution rates (Figure 5B). 
Decreased uptake rates due to low dilution rates resulted in 
activation of the PEP-glyoxylate cycle, for both growth on 
glucose and galactose. As galactose uptake is normally 
repressed in the wild type, we used the NagC mutant to enable 
derepressed uptake of galactose (Soupene et al 2003; Samir El 
et al 2009). At a dilution rate of 0.12 h" 1 , resulting in an 
uptake of 1.38 mmol gCDW -1 h -1 (g cell dry weight, gCDW), 
the NagC mutant exhibited a similar flux phenotype as was 
obtained for the wild type. In contrast to the wild type, 
however, the NagC mutant was capable to grow at a dilution 
rate of 0.35 h _1 on galactose, thereby achieving a galactose 
uptake rate of 3.53 mmolgCDW -1 h -1 (compared with a 
maximal rate of 2.00 mmol gCDW -1 h _1 for wild-type batch 
growth). At both increased galactose and glucose uptake rate, 
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Figure 5 Fraction of flux entering the PEP-glyoxylate cycle as a function of 
hexose uptake rate in batch (A) and chemostat (B) cultures. (A) Batch growth on 
galactose: the wild-type value is highlighted by a black circle. The PTS enzyme 
IIA glc mutant (crr) is highlighted in red. Values are the average of two independent 
experiments obtained by C-constrained flux analysis (Supplementary Tables 2 
and 3). For mutants with significantly altered flux distributions (labeled), a third 
independent experiment was conducted and the error bars represent standard 
deviations of these three independent experiments. (B) Chemostat growth on 
glucose and galactose: wild type on glucose (O), NagC mutant on galactose is 
highlighted in blue (♦); values from batch experiments are additionally listed 
( x ). Nomenclature: CS, chemostat culture; B, batch culture. Two biological 
replicates were conducted for each condition and error bars represent the 
deviation between the two experiments. 



the fraction of carbon entering the PEP-glyoxylate cycle flux at 
the acetyl- Co enzyme A branch point was significantly 
reduced. Thus, independent of the substrate, the fraction of 
acetyl-Coenzyme A that enters the PEP-glyoxylate cycle 
decreases with gradually increasing hexose uptake rate. As 
the PEP-glyoxylate cycle bypasses the NADPH-forming iso- 
citrate dehydrogenase of the TCA cycle, it adds potential 
metabolic flexibility to redox metabolism. Apparently, as the 
cell enters a condition where higher formation than consump- 
tion of NADPH occurs, e.g., during carbon-limited growth in a 
chemostat when relatively little biomass is synthesized, the 
PEP-glyoxylate cycle is activated, whereas the relative flux 
through the NADPH generating PP pathway remains unaltered 
(Figure 4A) . The PEP-glyoxylate cycle, however, is not the only 
mechanism to decouple catabolic carbon flow from NADPH 
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formation because E. coli also contains two actively used 
transhydrogenases that can interconvert NADPH and NADH 
(Sauer et al, 2004). 

Having established causality between uptake rate and the 
PEP-glyoxylate cycle, we wondered whether the decreased 
PEP-glyoxylate cycle observed in several mutants on galactose 
might be an indirect consequence of low cAMP levels during 
the increased catabolic flux regimes (Bettenbrock et al, 2007) . 
To test this hypothesis, we determined intracellular cAMP 
levels in wild type and mutants with a repressed PEP- 
glyoxylate cycle in galactose batch cultures (Figure 6) . In all 
mutants with increased growth rate on galactose, intracellular 
cAMP levels were significantly reduced compared with the 
wild type, reaching comparable reduced cAMP levels as 
measured for wild-type growth on glucose. The reduced PEP- 
glyoxylate cycle in Cra, IHF A, IHF B, Mlc and NagC is hence 
not a direct transcriptional control, but rather the indirect 
consequence of lowered intracellular cAMP levels that result in 
reduced Crp activity. Of the 36 transcription factors with target 
genes in TCA cycle, glyoxylate cycle and acetate secretion that 
could therefore potentially control the flux distribution at the 
acetyl-CoA branch point, thus only Crp actively controls flux 
during growth on galactose. 

Uptake of galactose is repressed in the wild type 

Somewhat surprisingly, several mutants grew faster and more 
than doubled the normal-specific galactose uptake rate of their 
E. coli wild-type parent (Figure 5A). Hence, we wondered 
whether E. coli might feature a similar non-optimal growth 
phenotype during exponential batch growth as was described 
for B. subtilis (Fischer and Sauer, 2005). As a measure of 
optimal growth, we plotted growth rate as a function of the 
experimentally determined biomass yield for the 91 mutants 
during growth on glucose and galactose. The product of both 
values represent the biomass productivity (gceiisghexoseh -1 ), 
and the thick line marks the biomass productivity that is 
achieved by the wild type (Figure 7) . 
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Figure 6 Intracellular cAMP concentrations relative to the wild type on 
galactose. Intracellular cAMP levels for mutants grown on batch galactose and 
the wild type grown on batch glucose relative to intracellular cAMP levels for the 
wild type grown on galactose. Intracellular cAMP levels were determined from 
two biological replicates for each strain from two separate shake flask 
experiments. Error bars represent the deviation between the two experiments. 



In contrast to B. subtilis, none of the E. coli mutants 
exhibited increased biomass productivity on its preferred 
substrate glucose, with the sole possible exception of Fis 
(Figure 7A), whose deletion leads to an increase in stable 
rRNA synthesis (Appleman et al, 1998). On galactose, 
however, the previously identified mutants with higher 
galactose uptakes (Figure 5A) clearly exhibited significantly 
increased biomass productivity, in particular Cra, IHF A, IHF B 
and NagC (Figure 7B), further supporting the conclusion 
that E. coli actively represses its galactose uptake rate at 
the expense of otherwise possible rapid growth. As also the 
PTS glc mutant Crr exhibited increased biomass productivity, 
repression of growth is directly or indirectly conducted by 
the PTS glc . 
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Figure 7 Growth optimality on glucose (A) and galactose (■). Growth 
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compared with the wild type (highlighted by a black circle) on (A) glucose and (B) 
galactose. Lines indicate equal biomass productivity (g (ghexoseh)" 1 ). Error 
bars represent standard deviations from at least three independent experiments. 
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Discussion 

From large-scale, 1 3 C- determined intracellular fluxes in 91 
transcriptional regulator mutants of E. coli, we identified active 
transcriptional control of hexose metabolism during two 
surprisingly distinct metabolic states. During maximum 
exponential growth, glucose metabolism was characterized 
through high metabolic rates with about fourfold higher 
overflow metabolism than respiratory TCA cycle fluxes. 
Galactose metabolism, on the other hand, was generally much 
slower, without overflow metabolism and with the PEP- 
glyoxylate cycle (Fischer and Sauer, 2003b) replacing the 
TCA cycle for respiration. Approximately 2/3 of the investi- 
gated regulators controlled absolute fluxes directly or indir- 
ectly through growth-related effects. Transcriptional control 
of flux partitioning, in sharp contrast, was confined to the 
acetyl-CoA branch point, where hexose catabolism can 
proceed via acetate overflow, the respiratory TCA cycle and/ 
or the respiratory PEP-glyoxlate cycle. The relative flux into 
the PP pathway, e.g., remained stable at ~30% over an order 
of magnitude of hexose uptake rates at very different biomass 
yields, demonstrating its primary catabolic rather than 
anabolic function in E. coli. 

On glucose, transcriptional regulation at the acetyl-CoA 
branch point controlled flux partitioning between the respira- 
tory TCA cycle and acetate overflow. A negative control of the 
TCA cycle/respiration branch was exerted by transcription 
factors ArcA, as was previously described (Perrenoud and 
Sauer, 2005), PdhR, and IHF A and B. All four are known to 
repress target genes in the TCA cycle and/or the respiratory 
chain (PdhR (Quail and Guest, 1995; Ogasawara et al, 2007), 
ArcA (Iuchi and Lin, 1988, 1992; Park et al, 1997; Liu and De 
Wulf, 2004; Perrenoud and Sauer, 2005; Vemuri et al, 2006), 
IHF A and B (Bongaerts etal, 1995; Green etal, 1997; Park etal, 
1997; Cunningham and Guest, 1998)). Increased capacity of 
the TCA cycle and/or respiratory chain in these mutants 
probably reduces acetate overflow. A positive control of the 
TCA cycle flux was exerted by the iron utilization regulator Fur, 
presumably through its known target genes in the TCA cycle 
(Escolar et al, 1998; Zhang et al, 2005). The potentially novel 
regulators GlpR, GlcC, QseB and HdfR could act either via 
positive control of TCA cycle/respiration or negative control of 
overflow metabolism. Our data do not allow drawing conclu- 
sions on their direct transcriptional control, because the flux 
alterations might also arise through indirect effects. 

On galactose, flux partitioning at the acetyl-CoA branch 
point was more complicated with the additional outgoing flux 
through the respiratory PEP-glyoxylate cycle and the TCA cycle 
flux mainly being confined to the anabolic purpose of 
supplying 2-oxoglutarate. Transcriptional control appeared to 
focus exclusively on the PEP-glyoxylate cycle, where we could 
distinguish two different groups of positively acting transcrip- 
tion factors. The only direct regulator was the cAMP receptor 
protein Crp (Kolb et al, 1993; Gosset et al, 2004), whose 
deletion entirely abolished PEP-glyoxylate cycle fluxes, as was 
described earlier for slow growing, glucose-limited chemostat 
cultures (Nanchen et al, 2008). Presumably, Crp induces 
expression of its target genes aceA and aceB (Gosset et al, 
2004; Zhang et al, 2005) upon allosteric activation 
through the high cAMP concentrations at low growth rates 
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(Bettenbrock et al, 2007). The second group consists of the 
indirectly acting transcription factors Cra, IHF A, IHF B, 
Mlc and NagC, whose deletion likewise abolished the 
PEP-glyoxylate cycle flux. Here the mode of action is the 
strongly increased galactose uptake rate and concomitantly 
decreased cAMP levels, which in turn deactivated Crp. Three 
lines of evidence support this hypothesis. First, we obtained 
similarly abolished PEP-glyoxylate cycle fluxes when increas- 
ing galactose uptake by deleting the glucose PTS-based 
repression of galactose uptake in the IIA glc mutant Crr 
(Meadow etal, 1990;Hogema etal, 1998). Second, modulation 
of galactose and glucose uptake rates through chemostat 
cultivations supports the inverse relationship between PEP- 
glyoxylate cycle flux and hexose uptake. This finding is 
consistent with the reported increased mRNA and protein 
levels of PEP carboxykinase and glyoxylate shunt enzymes at 
decreased glucose uptake rates (Ishii et al, 2007). Third, the 
low intracellular cAMP levels in our mutants with increased 
galactose uptake rates are fully consistent with the known 
inverse relationship between intracellular cAMP levels and 
catabolic rates (Matin and Matin, 1982; Notley-McRobb et al, 
1997; Ferenci, 2001; Fischer and Sauer, 2003b) and establish 
the link to Crp. Hence, we conclude that the only direct 
transcriptional regulator of the PEP-glyoxylate cycle on 
galactose is the cAMP-activated Crp protein, while the other 
five regulators act indirectly via cAMP and Crp. 

By selecting a large-set of transcriptional regulators, includ- 
ing all transcription factors known to regulate enzymes of 
central metabolism, we demonstrate the extent to which this 
flux splitting is subjected to transcriptional regulation. Overall, 
our absolute flux data demonstrate that control of flux splitting 
during growth on hexoses was confined to the acetyl-CoA 
branch point in E. coli. Of the 36 transcription factors known to 
target genes in pathways that diverge from the acetyl-CoA 
branch point, only one transcription factor on galactose and 
five plus potentially four others on glucose showed altered flux 
splitting. The primary focus of steady-state transcriptional 
control on the acetyl-CoA branch point, and thus the metabolic 
decision between the energetically efficient respiration and the 
less efficient but more rapid fermentation, was recently also 
demonstrated with only relative flux data for Saccharomyces 
cerevisiae (Fendt et al, 2010). Likewise, the identified networks 
of active transcriptional control of the glucose flux distribution 
were of similar size in yeast and E. coli with four and five (plus 
four) transcriptional factors, respectively. On galactose, how- 
ever, the active yeast network was quite extensive with more 
than 10 factors, while fluxes in E. coli where controlled by only 
one transcription factor. Why the yeast active transcription 
network of respiratory metabolism is much more extensive on 
galactose than the E. coli network remains unclear as a similar 
number of transcription factors target the TCA cycle (35 for 
S. cerevisiae (Monteiro et al, 2008) and 36 for E. coli (Keseler 
et al, 2005; Salgado et al, 2006)). It is tempting to speculate 
though that it is related to the fact that yeast does not invoke a 
novel respiratory pathway such as E. coli but rather needs to 
activate its TCA cycle (Fendt et al, 2010) . As different metabolic 
pathways were involved in E. coli galactose catabolism when 
compared with glucose, it was necessary to determine 
the much more quantitative and tedious absolute in vivo 
fluxes, which then enabled us to identify the single directly 
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controlling transcription factor Crp from the larger set of 
possible candidates. 

Previously, suboptimal biomass productivity was described 
for glucose-grown B. subtilis (Fischer and Sauer, 2003b). 
For batch glucose-grown E. coli, we demonstrate here that 
none of the investigated transcription factor mutants exhibited 
improved biomass productivity. On galactose in contrast, 
optimal growth appears to be transcriptionally repressed, 
because the Cra, IHF A, IHF B and NagC mutants grew much 
faster at almost unaltered biomass yields. As we could 
engineer a similarly improved phenotype by directly increas- 
ing the galactose uptake through removal of the glucose 
PTS-based repression with a Crr mutant, we provide evidence 
that E. coli actively represses its galactose uptake at the 
expense of otherwise possible rapid growth. 



Molecular Devices, Sunnyvale, CA). Specific growth rates were 
determined from log-linear regression of time-dependent changes in 
optical density from at least four data points during the exponential 
growth phase. 

Glucose, galactose and acetate concentrations were quantified 
enzymatically using commercial enzyme kits (Enzytec, Switzerland) . 
The uptake and secretion rates were determined from at least three 
technical replicates of two independent shake flask or microtiterplate 
experiments from two points (beginning of exponential growth and 
mid-exponential growth). The glucose, galactose and acetate concen- 
trations during exponential growth were plotted against the corre- 
sponding cell dry weights. Cell dry weight was calculated using a 
determined conversion factor of 0.41 and 0.51 gCDW) per liter per 
OD 6 oo for glucose and galactose, respectively. A linear fit was applied 
to calculate the slope. The inverse of the slope of glucose/galactose 
concentration against cell dry weight is the biomass yield in g mmol" 1 . 
The non-inversed slopes were further multiplied with the growth rate 
to get uptake and secretion rates. 



Materials and methods 

Strain and growth conditions 

All experiments were performed with E. coli BW25113 and its 
otherwise isogenic mutants from the Keio knockout library (Baba 
et al, 2006). Experiments were conducted with both independently 
generated mutants from the library. For clarity, the here used mutant 
nomenclature reflects the deleted genes (Supplementary Table 1). 
Batch physiology and flux experiments were performed in 35 ml 
cultures in 500 ml shake flask or 1.2 ml cultures in 2 ml 96-deep-well 
plates (Duetz et al, 2000) at 37°C, 300 r.p.m. and a shaking diameter of 
5 cm. Chemostat experiments were conducted in eight parallel 10 ml 
bioreactors as previously described (Nanchen et al, 2006). Briefly, 
aeration was achieved by pumping water-saturated air through the 
reactors and a constant temperature was assured through incubation 
in a 37°C water bath. A medium flow of 1 .5 ml h -1 and reactor volumes 
of 3.0, 4.3, 7.5 and 12.5 ml were used to obtain dilution rates of 0.5, 
0.35, 0.2 and 0.12 h" 1 , respectively. 

Frozen glycerol stocks were used to inoculate Luria-Bertani (LB) 
complex medium, supplemented with 50mgl _1 kanamycin for all 
mutants. All further cultivations were performed without antibiotics. 
After 6 h of incubation at 37°C and constant shaking, LB cultures were 
used to inoculate M9 minimal medium precultures upon overnight 
cultivation. The following day 0.3 ml and 10 |il of the M9 precultures 
were used to inoculate batch physiological and flux experiments for 
shake flask and 2 ml 96 deep-well plates, respectively. For chemostat 
experiments, the M9 precultures were used to inoculate 10-ml-scale 
bioreactors (Nanchen et al, 2006) with M9 minimal medium contain- 
ing 13 C-labeled substrate at a starting OD 600 of 0.05. The medium feed 
for 13 C-labeled glucose- or galactose-limited chemostat operation was 
initiated after 4-8 h of batch growth depending on the achieved 
biomass. 

The M9 medium contained, per liter of deionized water, the 
following: 0.8 g (NH 4 ) 2 S0 4 , 0.5 g NaCl, 7.5 g Na 2 HP0 4 2H 2 0 and 3.0 g 
KH 2 P0 4 . The following components were sterilized separately and 
then added (per liter of final medium) : 1 ml of 1 MMgS0 4 , 1 ml of 0.1 M 
CaCl 2 , 0.6 ml 0.1 M FeCl 3 6H 2 0, 2 ml of 1 mM filter-sterilized thiamine 
HQ and 10 ml of a trace element solution containing (per liter) 0.18 g 
ZnS0 4 7H 2 0, 0.12 g CuCl 2 2H 2 0, 0.12 g MnS0 4 H 2 0 and 0.18 g CoCl 2 
6H 2 0. Filter-sterilized glucose or galactose was added to a final 
concentration of 3gl _1 and lgl -1 for batch and chemostat experi- 
ments, respectively. For 13 C-labeling experiments, glucose or galactose 
was added entirely as the [1- 13 C] -labeled isotope isomere (>99%; 
Cambridge Isotope Laboratories, Andover, MA) or as a mixture of 20% 
(wt/wt) [U- 13 C] -labeled (>99%; Cambridge Isotope Laboratories) 
and 80% (wt/wt) of natural glucose or galactose. 

Physiological parameters 

Cell growth was monitored by determining the optical density 
at 600 nm (OD 60 o) using a spectrophotometer (Spectra Max Plus, 



Metabolic flux ratio analysis by GC-MS 

Aliquots of fractionally 13 C-labeled biomass were prepared as 
described previously for gas chromatography mass spectrometry 
(GC-MS) analysis (Fischer and Sauer, 2003a; Fischer et al, 2004; 
Zamboni et al, 2009). Briefly, cell pellets from 1 ml culture aliquots at 
an OD 600 of 0.7-1 .5 were collected by centrifugation and hydro lyzed in 
6M HC1 at 105°C for 24 h in sealed Eppendorf microtubes or 96-well 
PCR racks for shake flask and deep- well plate cultures, respectively. 
The hydro lysates were dried under a stream of air at ~60°C and 
then derivatized at 85°C in 20 (il dimethylformamide (Fluka, 
Switzerland) and 20 ul N-(tert-butyldimethylsilyl)-N-methyl-trifluoro- 
acetamide with 1% (vol/vol) tert-butyldimethylchlorosilane 
(Fluka, Switzerland) for 60 min (Fischer et al, 2004) . 

Derivatized amino acids were analyzed on a 6890 GC system 
(Agilent Technologies, Santa Clara, USA) combined with a 5973 Inert 
SL MS system (Agilent Technologies, Santa Clara, USA). The GC-MS- 
derived mass isotope partitioning of proteinogenic amino acids were 
then corrected for naturally occurring isotopes (van Winden et al, 
2002). The following eight metabolic flux ratios were calculated, as 
described previously (Fischer and Sauer, 2003a; Zamboni et al, 2005, 
2009; Nanchen et al, 2007): serine derived through glycolysis, serine 
derived through the ED pathway, oxaloacetate originating from PEP, 
PEP originating from oxaloacetate, oxaloacetate originating from 
the glyoxylate shunt, pyruvate originating from malate (upper and 
lower boundaries) and an upper bound on PEP derived through the 
PP pathway. 



Estimation of absolute fluxes by the 
13 C-constrained method 

Intracellular absolute carbon fluxes were estimated using the software 
FiatFlux with the previously described (Fischer et al, 2004; Zamboni 
et al, 2005, 2009) stoichiometric model that included all major 
pathways of E. coli central carbon metabolism, including the 
glyoxylate shunt and ED pathway. The reaction matrix consisted of 
25 unknown fluxes and 21 metabolite balances (including the 
three experimentally determined rates of glucose/galactose uptake, 
acetate secretion and biomass production). To solve the under- 
determined system of equations, the eight above metabolic flux 
ratios were used as additional constraints, as described before (Fischer 
et al, 2004; Zamboni et al, 2009). The first five ratios were used as 
equality constraints, while the latter three were used as boundary 
constraints. 

Fluxes into biomass were calculated from the known metabolite 
requirements for macromolecular compounds and the growth-rate- 
dependent RNA and protein contents (Emmerling et al, 2002) . Using 
the FiatFlux software (Zamboni et al, 2005), the sum of the weighted 
square residuals of the constraints from both metabolite balances and 
flux ratios was minimized using the MATLAB function fmincon, and 
the residuals were weighed by dividing through the experimental error 
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(Zamboni et al, 2009). All calculations were performed in Matlab 7.7.0 
(The Mathworks Inc., Natick, USA) . 



Estimation of absolute fluxes by the whole 
isotopologue modeling method 

High precision estimation of absolute fluxes was performed in selected 
cases by whole isotopologue balancing (Kleijn et al, 2005; van Winden 
et al, 2005). Briefly, we used cumomer balances and cumomer to 
isotopologue mapping matrices (Wiechert et al, 1999) to calculate the 
isotopologue partitioning of metabolites in a pre-defined stoichio- 
metric network model for a given flux set. The flux set that gives the 
best correspondence between the measured and simulated 3 C-label 
partitioning is determined by nonlinear optimization and denoted 
as the optimal flux-fit. All calculations were performed in Matlab 
7.7.0 (The Mathworks Inc.). The standard deviations for metabolic 
fluxes were derived from Monte Carlo simulations (Schmidt et al, 
1999). To mimic measurement errors, Gaussian noise was added 
to the measured 13 C-labeling data, and fluxes were re-estimated, 
thereby providing a measure for the sensitivity of the different 
fluxes. 



Intracellular cAMP measurements 

Cells were grown in batch cultures as described above. Of cultures 
growing in mid-exponential phase, 1 ml aliquots were taken in a 
constant temperature room kept at 37°C and were vacuum-filtered on a 
0.45 um pore size nitrocellulose filter (Millipore). Samples were 
immediately washed with two volumes of fresh M9 medium at a 
temperature of 37°C, containing the respective carbon source, and 
adjusted to pH of the culture at the time of sampling. After washing, the 
filter was directly transferred for extraction into 4 ml of 60% (vol/vol) 
EtOH/H 2 0 kept at 78°C for 2 min. After extraction, the solution was 
separated from the filters and centrifuged to remove residual cell 
debris and nitrocellulose. Cell extracts were thawed, dried at 120 |i bar 
and resuspended in 100 ul deionized H 2 0, of which 15ul were 
transferred into rubber-sealed HPLC tubes. Intracellular cAMP 
concentrations were determined by using an ion-pairing ultrahigh 
performance liquid chromatography-tandem mass spectrometry 
method (Buescher et al, 2010). Compounds were separated using a 
Waters Acquity UPLC with a Waters Acquity T3 end-capped reverse- 
phase column (150 x 2.1 mm x 1.8 um; Waters Corporation, Milford, 
MA, USA) and subsequently detected on a tandem mass spectrometer 
(Thermo TSQ Quantum Triple Quadropole with Electron-Spray 
Ionization, Thermo Scientific, Waltham, MA, USA) . 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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